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the wavefunctions and the spectrum are known to show fractal structures. Such 
wavefunctions are called critical. We present a phase diagram of the energy spectrum 
for varying the nonlinearity. It consists of three portions, a forbidden region, the 
spectrum of critical states, and the spectrum of stationary solitons. We show that the 
energy spectrum of critical states remains intact irrespective of the nonlinearity in the 
sea of a large number of stationary solitons. 
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1. Introduction 

The realization of Bose-Einstein condensation (BEC) in optical lattices has opened a 
new avenue for studying a variety of phenomena in condensed matter systems [BIS]- A 
major advantage of using ultracold atomic gases is that one can control the interatomic 
interactions and the lattice parameters in an extremely clean environment. This high 
degree of tunability enables one to study BEC in artificially designed structures which 
cannot be achieved in conventional solids. For instance, the one- dimensional bichromatic 
potential has been realized in a system of *^Rb atoms [3], and such a quasiperiodic 
potential has been studied theoretically [1] in the context of BEC. More recently, a new 
method for creating potentials through a holographic mask was introduced [5]. Using 
this technique, it seems feasible to experimentally generate exotic structures such as the 
Fibonacci lattice P-[8] and the Penrose tiling [9l[Tn] which are of interest as one- and 
two-dimensional quasicrystals, respectively. 

These low- dimensional quasicrystals have attracted considerable theoretical 
attention since there appear fractal wavefunctions which are neither extended nor 
localized [Hllll]. They are called critical states. The rapid progress in the study of BEC 
can set a stage for exploring these exotic states in experiments. The new ingredient that 
appears in the system of BEC is the nonlinearity caused by interatomic interactions 
which can be tuned by the Feshbach resonance. One might expect that fractal 
wavefunctions are fragile and are easily destroyed by the nonlinearity. Surprisingly, 
this is not the case. In fact, we demonstrate that there indeed exist critical states 
on the Fibonacci optical lattice. This is intended to stimulate experimental efforts to 
observe critical states in a cold-atom setuplil 

In order to describe a Bose-Einstein condensate on the Fibonacci optical lattice, 
we resort to the nonlinear Schrodinger equation with on-site potentials arranged in the 
Fibonacci sequence |13|. In the absence of the nonlinear term, it is known that all the 
eigenstates are critical, and that the spectrum shows a fractal structure |HlllI]- More 
precisely, the spectrum is singular continuous [TT1IT4] and called the Cantor spectrum. In 
order to elucidate the effect of nonlinearity on the critical states, we numerically solve 
the stationary nonlinear Schrodinger equation. As mentioned above, our numerical 
results show that the critical states persist despite the presence of the nonlinearity in 
the sea of stationary solitons p^ [T5]. With the aid of mathematical tools, we show that 
for any critical state, the "eigenenergy" must be included in the Cantor spectrum of the 
Schrodinger equation without nonlinearity. Further, we determine a forbidden region for 
the eigenenergy and the strength of the nonlinearity. Putting these together, we present 
a phase diagram of the energy spectrum for varying the nonlinearity. (See Fig. H] in 
Sec. O) The energy spectrum of the critical states retains its profile irrespective of the 
nonlinearity, while the number of the stationary solitons increases enormously as the 
nonlinearity increases. One might think that the presence of the sea of solitons makes 

X Quite recently, disorder effects were studied for transport in photonic quasicrystals [12]. Their 
experimental results show that certain disorder enhances the transport. 
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it difficult to experimentally detect the fractal profiles of the critical states. However, in 
the neighborhood of the forbidden region, such an experimental detection is expected 
to be possible. (See Sec. Eland Fig. Hlfor details.) 

Throughout the present paper, we will not treat dynamical properties of the 
wavefunctions for the nonlinear Schrodinger equation. However, we think we should at 
least stress that knowledge of the stationary states is not sufficient for understanding the 
dynamics such as diffusion of wave packets in quasiperiodic or random environment [TBI - 
[21] . For the random nonlinear Schrodinger equation, see a recent review [22] . 

The present paper is organized as follows: The precise definition of the nonlinear 
Schrodinger equation which we consider in the present paper is given in Sec. El In Sec. [31 
the numerical solutions of the model are obtained. In Sec. [H we apply the multifractal 
analysis to the solutions so obtained. The mathematical analysis for the model is given 
in Sec. [5l Our results are summarized as the phase diagram in Sec. [6l Section [3 is 
devoted to summary and conclusion. In Appendix A, we discuss effects of nonlinearity 
on localization. 



2. Preliminary 

Let us consider the one-dimensional chain which is determined by the Fibonacci rule. 
The £-th chain, S*^, £= 1, 2, . . ., consists of two symbols, A and B, and is constructed by 
the recursion, S^j^i = S^S^-i, with the initial condition, So = B and 5*1 = ^4: 

Si = A, S2 = AB, S3 = ABA, Si = ABAAB, ... 

We denote by A^ = the number of symbols in Se. Clearly, Ni is equal to the Fibonacci 
number because they satisfy N£+i = N^ + Ni^i with A'"o = A''i = l. In the limit j-oo, the 
Fibonacci chain is neither random nor periodic. In fact, it is quasiperiodic. 

The nonlinear Schrodinger equation for stationary states, {ipi,ip2, ■ ■ ■ ,iPn), on the 
Fibonacci chain is given by 

^(V'i+l+^^_l) + V^^^^ + ^|^^|'^i = ^^i for i = l,2,...,N, (1) 

where the hopping integral t and the coupling constant g are real, E is an 
"eigenenergy" [231425] . and the on-site potential Vi is given by 




Va if ^-th symbol of Se is A, 
Vb if i-th symbol of Se is B 



with real Va and Vb- We impose the Dirichlet boundary conditions, ipo = ipN+i = 0, 
or the periodic boundary condition, i/jn-h = '^Pi- We choose the normalization of the 
wavefunctions as J2iLi l'^iP = l- 
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3. Numerical analysis 

Using the shooting method, we numerically solve ([I]) with the Dirichlet boundary 
condition. We choose = — 1, Vb = +1 and t = l. For fixed E and g, we continuously 
vary the amplitude ipi so as to hit ipN+i = at site + A solution so obtained does 
not satisfy Yl!i=i iV^iP = 1 in general. So we set ipi = Xip'i with A > so as to satisfy 
YliLi iV'iP = 1- Then the coupling constant is given by g' = X^g. In the following we 
drop primes. 



0.2 - 








Figure 1. Numerical solutions of ([T]) are plotted in the E-g plane near the band center 
for chain length Ne — 377 with £—13. 

The numerical results for eigenenergies E are plotted in the E-g plane in Fig. [1] 
Each eigenenergy E is an increasing function of g. This is a consequence of the 
positivity of the nonlinear term. Unlike the linear Schrodinger equation, there appear 
many solutions whose number is beyond the total number of the sites of the chain 
N. These contain a large number of localized states. Since the localized states are 
due to the nonlinearity, we call them solitons. This phenomenon is already known as 
the appearance of localized states for lattice systems p^lTS] . There are solutions with 
different spatial profiles as shown in Figs. [21 It seems likely that all the states, extended, 
critical, soliton, and surface states, mix like a complicated "soup". However, we will 
show below that extended states are absent, whereas critical states and solitons exist. 
The surface states exist only for a system with the Dirichlet boundary condition. 

For the Harper model [26|[27] at criticality, effects of nonlinearity were discussed 
in |28]. Their numerical results show that the eigenmodes can be classified into two 
families: the conventional modes which are continuously connected to those in the 
corresponding linear model and new modes that have no origin in the linear model. If 
the latter modes are interpreted as stationary solitons due to the nonlinearity as in our 
case, their results of [28] are consistent with ours. 
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Figure 2. The left column shows critical wavefunctions in (a) a linear case and (b) a 
nonlinear case. The right column shows localized waveftmctions of (c) a surface state 
and (d) a soliton in a nonlinear case. All the data are for chain length Ni — 377 with 
£ = 13. 



4. Multifractal analysis 

As is well known, the scaling analysis which is called multifractal analysis [29|[30] is a 
very useful tool to determine whether or not a given wavefunction is critical [31] . Relying 
on this analysis, we numerically check the existence of critical states. We expect that 
the wavefunction which is shown in Fig. |2] (b) with 7^ leads to a critical state in 
the infinite-length limit of the chain. In the multifractal analysis, a wavefunction is 
characterized by /(a). For a wavefunction if) = {ipi, . . . jI/jm), the number of sites i 
satisfying \ipi\'^ ~ N~°' is assumed to be proportional to A^-/^«(°). The subscript £ and 
the chain length N are related through N = Ne with the Fibonacci number Ni whose 
definition is given in Sec. |2j We denote by /(a) the infinite-length limit of fe,{a). Using 
/(a), one can determine the character of a given wavefunction as follows: 

• An extended state shows a single point, (/, a) = (1, 1). 

• A localized state shows two points, (/,«) = (0,0) and (l,oo). 

• A critical state shows a sequence of smooth curves fe{a) which do not fall into the 
above two categories in the infinite- length limit. 

As to critical states, the corresponding sequence does not necessarily converge to some 
f{ct). But, it is known that /^(a) for a wavefunction with a multifractal character 
converges to a single smooth curve f{a). Our numerical results of fe{a) are shown in 
Fig. [3l The results show that the wavefunctions are critical because the behavior of 
the sequence of /^(a) is totally different from those of /^(a) for extended and localized 
states [31]. In consequence, the critical states persist in spite of the nonlinear term. 

However, it is still unclear whether or not the state shows perfect multifractality [32] 
within our numerical analysis since the right side of the profile of fe{a) is oscillating 
as i increases, and does not seem to converge to a single smooth curve in the limit 
£ — )■ cxD. The oscillation in Fig. [3] can be explained as the effect of the Dirichlet boundary 
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condition as follows: We treat only the case of the linear Schrodinger equation, i.e., ([T]) 
with g = We impose the periodic boundary condition. Let i/j = (t/'i, . . . ,iPn) and 
xj)' = {ip'i, . . . ,tp'^) be two independent eigenvectors with the eigenvalues, E and 
respectively. We assume that E is nearly equal to and that \ipi\^ ~ A^""^ and 
~ A^"°'i with ai > a[ > 0. Consider 

^ = ^- (2) 

Then, (f satisfies the Dirichlet boundary condition, (pi = (pN+i = 0. We assume that 
both of V' and i/^' show the same /^(a). In order to distinguish fila) of the wavefunction 
cp from that of if) or if)', we write fe{a) for fe{a) of cp. 

We want to show that fi^a) is not necessarily equal to fi{(y), and the deviation 
strongly depends on ai and a[. Consider the probability density at the site i. The 
contribution of the first term i/j in the right-hand side of ([2]) is written 

l^il' ~ iV""' (3) 
with some > 0. The contribution of the second term is written 



^1 ,/ 
^1 



with Aa := ai — a[, where we have assumed 

~ N-< (5) 

with some a'^ > 0. Therefore, we naively expect that the number of the site i satisfying 
~ A^~" is given by 

^/,(.) ^ |iV/^(") if /.(«)> /.(«-Aa), 
|^/,(a-Aa) if /,(«)< /,(a-Aa). 

However, the behavior in the second case is questionable. Clearly, from (jlj), the second 
case occurs for large a. Further, if a = + Aa > ai, then one has l^jjp ~ N~"\ In 
such a case, the second term in the right-hand side of ([2]) does not contribute to fi{a) 
although a- + Aa = a. From these observations, we conclude that iV-^^^""^") in the 
second case of ([6]) is reduced to some value N^'-^"'^ which satisfies 

In fact, as is well known, fi{a) must become a smooth curve having a single peak. To 
summarize, we obtain 

fi{a) = fe{a) if /^(a) > /^(a - Aa), 

Jiia) < fi{a) < fe{a - Aa) if fi{a) < fe{a - Aa). 

This explains the oscillation of fe{a) on the right side of the profile in Fig. [31 Thus the 
oscillation of the profile of the right side is due to the Dirichlet boundary condition. 
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Figure 3. (color online) fi{a) in (a) linear cases and (b) nonlinear cases. These 
numerical results show that the corresponding states are critical. Even in the linear 
case, ft{a) oscillates as ^ increases. 

5. Mathematical analysis 

We want to elucidate the spectral structure of Fig. [H and extrapolate it to the infinite- 
length limit of the chain (Fig. H]). 

First of all, it is practical to review the known results in the linear case, (7 = 0. 
Consider first the linear Schrodinger equation with the periodic boundary condition. In 
the infinite-length limit, all of the stationary states are critical and the energy spectrum 
becomes a singular continuous Cantor set. On the other hand, as to the Dirichlet 
boundary condition, the surface states appear and their eigenenergies form a pure point 
spectrum in addition to the Cantor spectrum of the critical states. 

In this section, we prove three theorems for the nonlinear Schrodinger equation: 
Theorem [1] states that there exists a forbidden region for two parameters, energy 
and nonlinearity, such that the nonlinear Schrodinger equation ([1]) with the periodic 
boundary condition has no solution. In the infinite-length limit, this forbidden region 
for the spectrum is common to both of the periodic and Dirichlet boundary conditions 
except for the spectrum of the surface states as we will show in Theorem El Theorem [2] 
states that an eigenenergy of a critical or extended state for the nonlinear Schrodinger 
equation is included in the spectrum in the case of (7 = in the infinite- length limit. 
This leads to the robustness of the critical states irrespective of the nonlinearity. 

Consider the linear Schrodinger equation, 

t{i)i+i + + Vii), = Eiji, (9) 

i.e., (II]) with g = and with the periodic boundary condition. We denote the Hamiltonian 
for (|9]) by Hq. We also denote the N eigenvalues by Sm, m = l,2, . . . , N, satisfying £j < 
for j < k. Define a set of real numbers, 

r := (-00, £1) u (^1 + ^, ^2) u • • ■ u {£n-i + g, £n) u {£n + g. +00), (10) 

for g > 0, where (■, ■) denotes an open interval. Here, ii g > £n\-i — then 
(^„, + (7, ^r^i) = for n = 1, 2, . . . , AT- 1. 
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Theorem 1. Let g > 0. Let E be an eigenenergy of the nonlinear Schrddinger 
equation ([1]) with the periodic boundary condition. Then, E^T. 

Proof. Let [ipi, . . . ^ip^) be a solution of the nonhnear Schrodinger equation ([1]) with 
the eigenenergy E. As in the assumption, we impose the periodic boundary condition. 
Then one can find n G {1, 2, . . . , A^} which satisfies £n< E < Snn, where £n+i '■= +00. 
Therefore, it is sufficient to show that E < £n + 9- Write Ui = and consider the 
linear Schrodinger equation, 

t{iJi+i + ^^Ji-l) + ViiJi + gUii)i = Ei)i, (11) 

with the additional potential gUi, and with the periodic boundary condition. Namely, 
we fix the additional potential gUi by using the solution {ipi, . . . ,iPn) of the nonlinear 
Schrodinger equation. Clearly, {ipi, . . . , ipjsf) is a particular solution of the equation ( fTTIl 
with the eigenvalue E. We denote by H the Hamiltonian for (fTTl) . and denote the N 
eigenvalues by £m, m = l,2, . . . , N, satisfying £j<£k for j < k. 

Note that the Hamiltonian H is written as H = Hg-\-gU with f/ = diag(f/i,f/2,. . . ,Un)- 
By the positivity of gll, we have Hq < if |§| As is well known, by applying the min-max 
principl^ to this type of an operator inequality, one can get an inequality between their 
eigenvalues. In the present case, we obtain f^r^f 1 < £^n+i • Combining this, the assumption 
E<Sn+.i and the fact that E is the eigenvalue of H, we obtain E<Sn- 

On the other hand, we have the bound, H <Ho+g. Applying the min-max principle 
again, we obtain Sn'^Sn + g- Combining this with the above result E<Sn, we obtain 
the desired bound E<£n+g. □ 

Thus the nonlinear Schrodinger equation ([1]) has no solution in the region 

r:={{E,g)\ EeT,g>0} (12) 

for two parameters, energy and nonlinearity. The region F is depicted as the non-colored 
region in Fig. HI by replacing "periodic" to "Dirichlet" in the proof. The forbidden 
region of the energy spectrum for the Dirichlet boundary condition will be treated in 
Theorem E] below. One might think that Theorem [T] is not too surprising because the 
deviation of the eigenenergy is less than or equal to g from t/j = iV^jp < 1. We stress that 
the statement of Theorem [1] includes that all the eigenenergies of the stationary solitons 
which are caused by the nonlinearity are also forbidden in the region F. The key idea 
of the proof is to introduce the potential gll into the linear Schrodinger equation. This 
enables us to apply the min-max principle to nonlinear problems for the first time. 

In order to determine whether a given state is critical or not, we must treat the 
infinite-length limit. Let {N{k)}'^^ be an increasing sequence of the length N = N{k) 
of the present chain. Let '0*^'^-' = ip2^\ . . . , -ijj^^j^f^-^) be a solution of the nonlinear 

§ When two self-adjoint operators, A and B, satisfy {(p^Aip) < {if,B(fi) for any state (p, we write 
A<B. 

II See, e.g., Sec. XIII. 1 in the book [33] , 
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Schrodinger equation ([T]) with the eigenenergy E^'^'^ for the chain with the length N{k). 
Write uf^'^ = IV'-'^^P, and [/max = niaxj f/f'^^ We assume that hni^^^ Umax exists. If not so, 
we take a subsequence. We say that the wavefunction 'ijj^''^ is locahzed if limk-^ Umlx > 0, 
and t/?*^*^) is critical or extended if linifc.^ ?7m2x = 0. We write if^^^^^ for Ho of ([9]) with 
the length N{k) and with the periodic boundary condition. The spectrum of -f^(fjv(fc) is 
given by a{H^^^^^^):={Si, . . . ,^iv(fc)}. 

Theorem 2. Consider the nonlinear Schrodinger equation with periodic /Dirichlet 
boundary condition. If a sequence {t/'*^*^^}^^ of the solutions is critical or extended in 
the infinite-length limit of the chain, k^oo, then the distance between the eigenenergy 
^(k) of ^l^^'^^ and the spectrum a^H^j^f^j^-j) must go to zero in the limit /c— J-oo. 

Proof. First consider the Hamiltonian H^^f^-^ for the linear Schrodinger equation ( ITTi) 
with the Dirichlet boundary condition in the proof of Theorem [1] with the additional 
potential U^'''' and with the length N{k). Clearly, '0'^'^) is the eigenvector of H^^^^^^ with 
the eigenvalue E^'^K This yields 



In order to obtain the lower bound for the right-hand side, we introduce the system of 

',N(k)- 



the complete orthonormal eigenvectors {u^^^^^} for H^j^,j^y Namely, 



with the eigenvalues Ei (k). Then, one has the expansion, 

^W^^^W^W ^ith ^|af)|' = l. (15) 

i i 

Using this, the right-hand side of f|T3|) is evaluated as 

>, (e(^) V'^'^)) (e(^) - 4'=))' 



(fc) 

a) 



(16) 



>m^in|(EW-4'=))'|. 
Substituting this into the right-hand side of (fT3|) . we obtain 

{h^w - H^Mk))' ^^'^) > [dist(a«^(,)),EW)]'. (17) 



Next, let us obtain an upper bound for the left-hand side of flTTl) . Note that 
HN(k)-Ho,N{k)=9U^''^+AHo^N(^k), where AHo^N(^k) = H^j^^^^-H^j^^^^^ i.e., the deference 
between the Hamiltonian Hq (Q with the Dirichlet and with the periodic boundary 
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conditions. Clearly, we have 



(fc) 



Since ^ff/^^^^ and Aifo,Af(fc) are self-adjoint, Schwarz inequality yields 
Substituting this into the right-hand side of flTSj) . we have 



+ 2 



< 



■n 2 



+ |tK/l^f¥ + l<\)P 



Combining this with the above bound ( |T7l) . we obtain 



dist(a«^(,)),E(^))<|(7 



(19) 



(20) 



(21) 



If i/j^*^) is critical or extended, then this right-hand side is vanishing as k-^oo. 

Clearly, the same statement holds with the periodic boundary condition. In 
this case, the second term in the right-hand side of (I2T1) does not appear since 



AH, 



0,N{k) '■ 



-0. 



□ 



We recall the well known fact that the spectrum (j{Hq^) of the linear model in 
the infinite- length limit is singular continuous and has zero Lebesgue measure [6| [TTllTi] . 
Theorem [2] states that all of the eigenenergies of critical or extended states in the 
nonlinear model fall into the set a{HQ^). This implies that the spectrum of critical or 
extended states in the nonlinear model in the infinite-length limit is a subset of (7{Hq^), 
and has zero Lebesgue measure. In general linear models, a set of extended states is 
defined to have a spectrum having nonvanishing Lebesgue measure. Therefore, in this 
sense, there is no extended state in the present nonlinear model. However, we cannot 
conclude, from Theorem |2l that critical states indeed exist in the nonlinear model, and 
that all of the critical states in the linear model survive switching on the nonlinearity. 
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This expectation is supported by our numerical results. Actually, as shown in Fig. [H 
each of the critical states is continuously connected to that in the linear model for varying 
the strength of the nonlinearity. We also remark that eigenenergies of surface states due 
to the Dirichlet boundary condition can appear outside the spectrum cr(if^^^^^) in 
general. 

The quantity f/j = \'4'i\^ which we introduced in the proof of Theorem [1] can be 
interpreted as an effective potential due to the nonlinearity. Since Ui is vanishing for 
critical or extended states in the infinite-length limit, one might think that Theorem [2] 
is a trivial consequence of this fact. This is not true because the effect of the 
nonlinearity for the whole chain is estimated by g Ui. Using the normalization 
condition = 1, one has gYli^i — 9^i\i^i\'^ = 9- Thus, for a fixed g, the 

effect of the nonlinearity is 0{1) irrespective of the length of the chain. We numerically 
check that the difference between critical wavefunctions il^g=o for g = Q and il^g^ for g^Q 
is 0{1) in the sense of norm | — '09=o| I irrespective of the chain length. Here, V^t*) 
is continuously connected with VcK)- Fo^' ^ typical V'st^); we have | |V'g7«i — V'cf=o| I —0.208 
for g ~ 0.098 and 1= 13, where I is defined in Sec. [2j Here the eigenenergy of ■i/'g=o is 
E~-0.187. 

Roughly speaking, nonlinearity does not change the character of eigenstates of a 
linear Schrodinger equation. Since the property of the on-site potential is not used in the 
proof of Theorem [21 one can expect that nonlinearity does not change the localization 
character of eigenstates of a random linear Schrodinger equation, too. Actually, we can 
justify this type of statement for a certain class of nonlinear Schrodinger equations. The 
precise statement of the theorem and its proof are given in Appendix A 

We write Too for the forbidden region F of ( IT2l) in the infinite-length limit, i.e., 
Foo := limAT^oor- Theorem |3] below states that the forbidden region for the spectrum 
of ([1]) with the Dirichlet boundary condition is identical to the forbidden region for 
the periodic boundary condition in the infinite-length limit except for the spectrum of 
the surface states. In other words, the spectrum of the surface states can appear on Foo- 

Theorem 3. Let g >0. Let {i/^^'^^'^i be a sequence of the solutions of the nonlinear 
Schrodinger equation ([1]) with the Dirichlet boundary condition and with the eigenenergy 
E^^^ such that 



lim 



(k) 



lim 

A;— ^-oo 



(k) 
N{k) 



0. 



Then, 



lim dist(M\F,E 



0, 



(22) 



(23) 



where F is given by f lTOj) andM.\r is the complement ofT. Namely, in the infinite- length 
limit of the chain, ([T]) has no solution in the forbidden region Too except for the surface 
states which are localized at the surface {1, A^}. 

Proof. Write H^^j^-^ for the Hamiltonian H in the proof of Theorem [1] with the additional 



potential uf^^ = and with the chain length N = N{k). Then the min-max principle 



',(fc)|2 
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yields a{H^^i^~^) nr = for any k. On the other hand, the same argument as in the proof 
of Theorem [2] yields 



dist(a(i7^(,)),E«)<|t|Vl^! 



',(fe)|2 



12 



(24) 



From the assumption (122|) . this right-hand side is vanishing as k^oo. These imply the 
desired result (1231) . □ 




E 



Figure 4. (color online) The spectral structure which consists of the Cantor spectrum 
of critical states (vertical red lines), the sea of stationary solitons (gray), and the 
forbidden region, Too (non-colored) , in E-q plane with the periodic boundary condition. 
An example of a region where experimentally detectable critical states exist are marked 
as stars. 



6. Phase diagram of the energy spectrum 

Let us describe the spectral properties of the nonlinear Schrodinger equation ([1]) in 
the limit £ — )■ oo. A schematic phase diagram is shown in Fig. HI It consists of three 
portions, the Cantor set for the critical states, the forbidden region T^, and the soliton 
sea. From the numerical results, we confirm that there exist critical states with the 
finite coupling g (Fig. [3] (b)). Combining this with Theorem [2l the Cantor set for the 
critical states is shown as vertical (red) lines in Fig. HI We numerically found solitons 
due to the nonlinear effects (Fig. [2] (d)). The point spectrum of the solitons must be 
distributed outside the forbidden region Too from Theorem [1] and |3l We cannot exclude 
the possibility that an eigenenergy of a soliton hes just on the Cantor set. It is obvious 
that surface states are absent with the periodic boundary condition, and not shown in 
Fig.H 



7. Summary and conclusion 

We studied the stationary states for the nonlinear Schrodinger equation on the Fibonacci 
optical lattice. We found that the nonlinearity does not destroy the critical states which 
exist in the absence of nonlinearity and exhibit fractal properties. The existence of these 
states was confirmed numerically using multifractal analysis. To our knowledge, this 
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kind of analysis is applied to the field of BEC for the first time. We also showed that 
the energy spectrum of the critical states remains intact irrespective of the strength 
of the nonlinearity. Besides the critical states, there is a large number of localized 
solutions, solitons, resulting from the nonlinearity. These solitons may seem an obstacle 
to observing critical states. However, we found the forbidden region for solitons, in the 
neighborhood of which the experimental detection of critical states is expected to be 
possible. Our analysis is intended to stimulate such an experimental effort to observe 
exotic critical states in optical lattices. 

The nonlinear Schrodinger equation ([1]) is nothing but the discrete Gross-Pitaevskii 
equation for BEC. Therefore, the chemical potential is equal to some eigenenergy E. 
In real experiments, the controllable parameters are the total number of the particles 
and the coupling constant g. In the nonlinear Schrodinger equation ([T]), these two 
parameters appear as a single parameter g which is the effective coupling constant under 
a normalization condition of the wavefunctions. Within a mean field approximation for 
many-body BEC systems, the effective single-body state which has the lowest internal 
energy is most likely to be realized as the ground state. We numerically checked, 
for relatively small effective couplings g (up to 1), that the ground state is given 
by the eigenstate of ([1]) with the lowest eigenenergy E. Then, E is identical to the 
corresponding the chemical potential /i. When a chemical potential /i is given instead of 
the total number of the particles, we can also numerically determine the internal energy. 

The method presented in this paper equally applies to a higher dimensional model 
in which the on-site potential in each direction is arranged by a generic quasiperiodic 
rule such as the Fibonacci rule. Actually, the eigenstates for the corresponding 
linear Schrodinger equation have a product form of the one-dimensional eigenstates. 
(See, for example, [331 - I36] .) This implies that, in order to detect critical states 
or fractal wavefunctions on BEC, experimentalists do not necessarily need to stick 
to a one- dimensional system. Also our method is applicable to a wide class of 
nonlinear Schrodinger equations for studying the properties of stationary states. For 
example, we can treat bichromatic on-site potentials |1] which is considered as the 
Harper equation |26[l27]. two different hopping integrals arranged in the Fibonacci 
sequence [8|l32]. and different types of nonlinearities such as the Ablowitz-Ladik one |37] . 
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Appendix A. Nonlinear effects for localization 

Although Theorem H] below holds for a wide class of nonlinear Schrodinger equations 
which have a localization regime in the spectrum of the corresponding linear Schrodinger 
equation, we consider a random nonlinear Schrodinger equation in one dimension as a 
concrete example. The nonlinear Schrodinger equation is given by replacing the on-site 
Fibonacci potential with a random potential in ([1]). As is well known, all the eigenstates 
in the corresponding linear system are localized for a general class of randomness in one 
dimension]^ As to the nonlinear eigenvalue problem, a certain set of localized eigenstates 
is proved to exist for more general setting [23ti25] . See also a related article [39] . 

In the following, we will prove that there is no stationary solution of the nonlinear 
Schrodinger equation such that the solution exhibits conventional properties of critical or 
extended states. Unfortunately, we cannot exclude the existence of certain pathological 
states which are not localized. We believe that such pathological states cannot appear 
for standard systems. Thus, the localization of the stationary states are expected to 
survive switching on the nonlinearity |30]. Clearly, our result is consistent with the 
previous results p3l425] . 

Let {A(fc)}^]^ be a sequence of finite lattices satisfying 

A(l) C A(2) C ■■■ C A(A;) C ■■■ . 

Let 'j/''''^^ = . . . be a stationary solution of the nonlinear Schrodinger 

equation on the lattice A(fc) with the eigenenergy E^^\ where N{k) is the number of the 
sites in k{k). We choose the sequence {k{k)}'^^i so that the eigenenergy E'^'^^ converges 
to some value E in the infinite- volume limit k ^ oc. Set 

ci>^'^^ := (0f ), := (A.l) 

maxj \%pi '\ 

For the wavefunction i/j^'^), we introduce the following two conditions: 

hm hm V Uf ) ' = oo, (A.2) 

fc'— >oo fc— >oo I 
iGA(fc') 



and there exists a finite lattice VL such that 



lim inf max 



> 0. (A.3) 



The former condition (\A.2\\ implies that the wavefunction i/?*^*^^ is not localized, and it 
does not split into two portions, a localized part and the rest, such that the distance 
between two portions becomes infinity in the infinite- volume limit. In fact, if 

lim lim Uf^ < oo, (A.4) 

fc'— >oo A;— >oo ' ^ I 
iGA(fc') 

^ See, e.g., the book [38] . 
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then 0^'^) converges to a localized state in the infinite- volume limit even for a critical or 
extended state -0*^^^ 

The latter condition ( lA.Sp implies that the wavefunction (p'^^'^ does not disappear 
from finite regions. 

We also consider the corresponding linear Schrodinger equation, and write Ho^N{k) 
for the Hamiltonian. The statement is given in a generic form as follows: 

Theorem 4. Let i/?'-^-' be a stationary solution of the nonlinear Schrodinger equation. 
Suppose that the eigenenergy E in the infinite-volume limit is an interior point of 
the localization regime in the energy spectrum of the corresponding linear Schrodinger 
equation. Then the solution xl)^^'^ cannot simultaneously satisfy the above two conditions 
M.^) and M.3j) . In other words, if a stationary solution is purely critical or extended 
in the sense of ^A.B^) and satisfies the energy condition, then its dominant part in the 
sense of the absolute value of the wavefunction cannot appear in any finite region. 



Proof. Assume that i/?^'^-' satisfies the conditions (]A.2p and (]A.3I) . From this assumption, 



we will show that one can construct an extended or a critical state for the corresponding 
linear Schrodinger equation in the infinite- volume limit. 

From (IA.3|) . there exist a subsequence {A(/cj)}°^^ of {A(A;)}^^ and a site G 



such that ^j-^^'' converges to some nonzero value ^j-^'' > as j — i- oo. Therefore, we 
can obtain a wavefunction (^(°°) in the infinite-volume limit by using the diagonal trick 
around the site Iq. Clearly, the wavefunction (j)^°°'> is nonvanishing, and non-normalizable 
from the condition flA.2p . 

In the same way as in ( ITTi) . we define by -ffAr(fc) the Hamiltonian with the additional 
potential U which is determined by ■ Then, we have ^^(^^■^(f)^'''^ = E^^^(j)^^K 
From these, we have 



J-KX) 



E(x(°°),0(°°)>, (A.5) 



where is a subsequence of {kj}JLi, and x^''^ is a function which converges to a 

rapidly decreasing function x^"^^ in the limit j— j-oo, and Hq oo is the Hamiltonian of the 
linear Schrodinger equation in the infinite- volume limit. This result implies that 0*^°°^ 
is a generalized eigenvectoiE] of Hq ^o- Since 0^°°^ is non-normalizable, this contradicts 
with the assumption that E is an interior point of the localization regime. □ 

The result ^B) implies that the vector = (..., 0S°°V ■ ■ , 4°°^ ■ ■ ■ ) 

locally satisfies the linear Schrodinger equation as 

t (^S+i^ + 0S) + V^^^ = Ecl>^\ (A.6) 

If there exists a localization regime in the corresponding linear Schrodineger 
equation, then the statement of Theorem H] generally holds. For example, the 

+ See, e.g., Sec. 4 of Chap. 1 in the book [41] . 
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Harper model has a localization regime for certain values of the coupling constant 
of the quasiperiodic potential. Therefore, there is no stationary solution which is a 
conventional critical or extended state by switching on nonlinearity under the condition 
that the nonlinear potential U vanishes in the infinite- volume limit. As a result, we can 
expect that the localization of the stationary states in the model is not destroyed by 
nonlinearity. 

Finally, we remark the following: Localization of stationary states may not lead 
to dynamical localization for nonlinear systems. (For related articles, see [T61 - I2T1H2] .) 
However, the existence [231425] of stationary localized states implies that, if nonlinearity 
dynamically destroys localization in a linear system, whether a wavefunction is 
dynamically localized or not strongly depends on the initial wavefunction. Actually, if 
an eigenstate is localized as a stationary solution of the nonlinear Schrodinger equation, 
then the time evolution of the state is localized, too. 
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